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Abstract 

The LVD detector, located in the INFN Gran Sasso National Laboratory (Italy), 
studies supernova neutrinos through the interactions with protons and carbon nu- 
clei in the liquid scintillator and interactions with the iron nuclei of the support 
structure. We investigate the effect of neutrino oscillations in the signal expected in 
the LVD detector. The MSW effect has been studied in detail for neutrinos travel- 
ling through the collapsing star and the Earth. We show that the expected number 
of events and their energy spectrum are sensitive to the oscillation parameters, in 
particular to the mass hierarchy and the value of ^13, presently unknown. Finally 
we discuss the astrophysical uncertainties, showing their importance and comparing 
it with the effect of neutrino oscillations on the expected signal. 
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1 Introduction 

There are many experimental works suggesting neutrino conversion among 
flavors in the recent few years, through the study of atmospheric [1], solar [2] 
[3] [4] [5] [6] [7], reactor [8] and accelerator [9] neutrinos. The interpretation of 
all these phenomena in terms of neutrino oscillations is rather robust, because 
it is able to include all the experimental data (except the "not yet confirmed" 
LSND [10] signal), even if the expected oscillatory behavior, in terms of the 
observable L/E, has not been yet experimentally observed (preliminary results 
that show a low significance hint for a oscillatory behavior have been found by 
a re-analysis of the SK data [1]). An interesting fact is that the inclusion of 
the MSW effect [11] permits a consistent interpretation of KamLAND results 
and the 'high energy' solar neutrino data [5,6,7]. 

In the standard three flavor scenario, six parameters must be determined by 
oscillation experiments: 3 mixing angles (^^soh ^13, ^^atm )i 2 squared mass dif- 
ferences (AmgQj and Am^^j^) and 1 CP- violation phase 5. A recent analysis 
of all the available experimental data [12] constrains the "atmospheric" and 
"solar" parameters to be in the following 99% C.L. ranges (compare also with 
the results in [13]): 



Oscillal ion i")arMin('i(\r 




ccutral valiu^ 
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solar mass splitting 




= (8.0 ± 0.3) 10-5 




(7.2 H- 8.9) lO-^eV^ 


atm. mass splitting 


Amatm 


= (2.5 ± 0.3) 10-^ 




(1.7-=- 3.3) 10-=^ eV^ 


solar mixing angle 


tan^ 6*801 


= 0.45 ± 0.05 




30° < ^soi < 38° 


atm. mixing angle 


sin^ 26'atm 


= 1.02 ± 0.04 




36° < ^atm < 54° 



However the other parameters are not completely determined: the ^^13 mix- 
ing angle is only upper limited, mainly by the Chooz experiment data [14] 
(sin^6'i3 < 3. 10~^ at the 99% C.L.), the sign of Am^t^, (that flxes the so- 
called mass hierarchy) is completely unknown, as well as the CP-violation 
phase 5. 

* Corresponding author: Marco Selvi, c/o INFN - Sezione di Bologna, via Irnerio 
46, 40126 Bologna - Italy, email: selvi@bo.infn.it 
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Because of the wide range of matter density in the stellar envelope, a super- 
nova explosion represents a unique scenario for further study of the neutrino 
oscillation mixing matrix. Indeed neutrinos can cross two resonance density 
layers and therefore the resulting possible mixing scenarios are different from 
the solar ones. The emerging neutrino spectra are sensitive to the sign of 
^^atm to the value of ^13. 

Before proceeding, it is important to recall that, at present, there is not a 
unique theory of supernova explosions. Till now, numerical investigations of 
the "standard model" based on a delayed scenario of the explosion [15] failed to 
reproduce the explosion. On top of that, other models are being studied where 
rotation [16] or magnetic field [17] play an essential role. In the following, we 
will use a simple description of the neutrino flux that does not contradict the 
SN1987A events seen by Kamiokande-II [18], 1MB [19] and Baksan [20], see 
e.g. [21] for a discussion, although it is not able to take into account the events 
seen in Mont-Blanc observatory [22,23]. This "standard" description, however, 
corresponds to the expected neutrino emission in the delayed scenario and in 
the last phase of the collapse with rotation [24]. For this reason, we take it as 
a useful starting point for the investigation of the impact of oscillations in the 
neutrino signal from a supernova. 

The main aim of this paper is to show how neutrino oscillations affect the 
signal detected by the LVD observatory in the INFN Gran Sasso National 
Laboratory, Italy. We also evaluate the impact on the signal of the astrophysi- 
cal parameters of the supernova explosion mechanism, such as the total energy 
emitted in neutrinos, the star distance, the neutrino-sphere temperatures and 
the partition of the energy among the neutrino flavors. 

In section 2 we describe the characteristics of the neutrino fluxes emitted dur- 
ing a gravitational core collapse. In section 3 the neutrino oscillation mecha- 
nism is shown, in particular the peculiarities of the MSW effect in the super- 
nova matter and in the Earth. The LVD detector and the relevant neutrino 
interactions both in the liquid scintillator and in the iron support structure 
are described in section 4. The impact of neutrino oscillations in the signal 
expected in the LVD detector is presented in section 5 while the uncertainties 
in the astrophysical parameters and their effect on the results are discussed 
in section 6. Finally, the conclusions are drawn in section 7. Two appendices 
complete this work, describing in more detail the MSW effect calculation in 
the Earth (A) and the neutrino interaction with the iron of the LVD support 
structure (B). 

Preliminary results have been presented previously in [25], [26] and [27]. 
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2 Supernova neutrino emission 



At the end of its burning phase a massive star (M > SMq) explodes into a 
supernova, originating a neutron star which rapidly cools down (order of tens 
of seconds) emitting about 99% of the liberated gravitational binding energy 
in neutrinos. 

The time-integrated spectra can be well approximated by the pinched Fermi- 
Dirac distribution. For the neutrinos of flavor a, we have 



where D is the distance to the supernova, E is the neutrino energy, 
is the time-integrated energy of the flavor I'a, represents the effective 
temperature of the i/a gas inside the star, r)a is the pinching parameter, 
F^^rja) = /o°° x^/(e^~^" + 1) dx is the normalization factor. In most of this 
work we assume for simplicity that = for all neutrino flavors; this choice 
results in the relation {Ea) — 3.15 between the mean neutrino energy 
temperature. 

Due to different trapping processes, the neutrino flavors originate in layers of 
the supernova with different temperatures. The electron (anti) neutrino flavor 
is kept in thermal equilibrium by P processes up to a certain radius usually 
referred to as the "neutrino-sphere", beyond which the neutrinos stream off 
freely. However, the practical absence of muons and taus in the supernova 
core implies that the other two neutrino flavors, here collectively denoted by 
{i^/i, i^ii, ^t), interact primarily by less efficient neutral-current processes. 
Therefore, their spectra are determined at deeper, i.e. hotter, regions. In ad- 
dition, since the content of neutrons is larger than that of protons, ly^s escape 
from more external regions than u^s. This rough picture leads to the hier- 
archy {El, J < {Ej^J < {Ey^). Typical ranges for the average energies of the 
time-integrated neutrino spectra obtained in simulations are {E^^) = 10 — 12 
MeV, {Ep^) = 11-17 MeV, and (E^J = 15-24 MeV [28,29]. However, 
recent studies with an improved treatment of v transport, micro-physics, the 
inclusion of the nucleon bremsstrahlung, and the energy transfer by recoils, 
find somewhat smaller differences between the Pg and spectra [30]. 

The amount of the total binding energy Ei, taken by each flavor is La = fu^ Eb, 
with ^ = 17 - 22%, /p^ = 17 - 28%, = 16 - 12% (see e.g. [31] ). Thus, 
the so-called "energy equipartition" has to be intended as "within a factor of 
two" [30]. 

In the following, if not specified differently, we assume a galactic supernova 
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Fig. 1. Neutrino energy spectra at the neutrino-sphere. 



explosion at a typical distance of D = 10 kpc, with a total binding energy 
of Eh = 3 ■ 10^^ erg and perfect energy equipartition f^^ = fp^ = f^^ = 1/6. 
We also assume that the fluxes of u^, Ur, i^n, and Pr are identical; we fix 
T^jTp^ = 1.5, T^^jTp^ = 0.8 and Tp^ = 5 MeV [30]. With these assumptions 
the resulting neutrino energy spectra generated inside the star are shown in 
figure 1. 



We postpone to section 6 a discussion about the implications to the expected 
number of events in the LVD detector due to the uncertainties in the astro- 
physical parameters. 
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3 Neutrino flavor transition in the star and in the Earth. 

In the study of supernova neutrinos, and Vt are indistinguishable, both 
in the star and in the detector, because of the corresponding charged lepton 
production threshold; consequently, in the frame of three-flavor oscillations, 
the relevant parameters are just (Arrigoi, UI2) and (Am^^j^^, t/^g) ^ . 

We will adopt the following numerical values: Amg^i = 8 • 10~^ eV^, Am^^^ = 
2.5 ■ 10~^ eV^, = 0.33; the selected solar parameters (Am^Qj, UI2) describe 
the LMA solution, as it results from a global analysis including solar, CHOOZ 
and KamLAND v data [12]. 

As described in figure 2 neutrinos, in the normal mass hierarchy (NH) scheme, 
cross two so-called Mikheyev-Smirnov-Wolfenstein [11] resonance layers in 
their path from the high density region where they are generated to the lower 
density one where they escape the star: one at higher density (H), which 
corresponds to (Am^^j^j, ^/g^g) and p = 300 -=- 6000 g/cm^ ^, and the other 
at lower density (L), corresponding to {Aml^i, UI2) and p = 5 100 g/cm^. 
Antineutrinos do not cross any MSW resonance [32,33,34]. 

For inverted mass hierarchy (IH), transitions at the higher density layer occur 
in the P sector, while at the lower density layer they occur in the v sector. 

Neutrinos arc originated in regions of the star where the density is very high, 
so that the effective mixing matrix in matter is practically diagonal. Thus the 
created neutrino fiavor eigenstate is completely projected into one neutrino 
mass eigenstate (represented by the thick purple line in figure 2). Then the 
neutrino starts its path through the matter to escape the star. If the matter 
density changes in a smooth way, then the propagation is said to be "adia- 
batic" . It means that the neutrino propagates through the star being the same 
mass eigenstate (i.e., referring to figure 2, staying over the same thick purple 
line). The adiabaticity condition depends both on the density variation and 
on the value of the oscillation parameters involved. 

Given the energy range of supernova v (up to ~ 100 MeV) and considering 
a star density profile p oc 1/r^, the L transition is adiabatic for any LMA 
solution values. Thus the probability to jump onto an adjacent mass eigenstate 
(hereafter called flip probability) is null [Pl = 0). The adiabaticity at the H 
resonance depends on the value of U^^ in the following way [34]: 

Ph oc exp [- const U^^ {AmlJEf'^ ] 

^ = cos^ • cos^ O12 ~ cos^ 6112, = cos^ • sln^ 9i2 ^ sin^ 612 and U^^ = 

^ the values are respectively for Ei, equal to 100 and 5 MeV 
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Fig. 2. Crossing level scheme for Normal (top) and Inverted (bottom) hierarchy. Solid 
thick purple lines show the eigenvalues of the effective Hamiltonian as function of 
the electron number density. The thin black lines correspond to the energy of the 
flavor levels and Vx- Negative values of the electron number density are related 
to the antineutrino channel. 

where Ph is the flip probability at the H resonance. 

When [/^g > 5 ■ 10"'' the conversion is completely adiabatic {ad) and the flip 
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probability is null (Ph = 0); conversely, when [/^^ < 5 ■ 10^^ the conversion 
is completely non adiabatic (na) and the flip probability is Ph = 1- We used 
in the calculation U^^ = 10~^, which is just behind the corner of the CHOOZ 
upper limit, for the adiabatic case and U^^ = 10~^ for the non adiabatic one. 

For neutrinos, in the NH-ad case z/g generated in the star arrive at Earth as 
z/3, so their probability to be detected as z/g is L^fg ~ 0. Thus, the detected 
come from higher-energy i/^ in the star that get the Earth as 1^2 and i^i. 
If the H transition is not adiabatic or if the hierarchy is inverted the original 
Ue get the Earth as U2 and their probabihty to be detected as Ue is U^2 ~ 0.3. 

For antineutrinos, in the NH case or in the IH-na, the z/g produced in the 
supernova core arrive at Earth as i^i, and they have a high {U^^ ~ 0.7) proba- 
bility to be detected as Pg- On the other hand, the original arrive at Earth 
as z/2 and z/g and arc detected as Pg with probability U^2- 
In the lH-ad case the detected Pg completely come from the original, higher- 
energy Pa; flux in the star. 

The oscillations scheme can be summarized as: 



for inverted hierarchy, 

where F°„^ are the original neutrino fluxes in the star and F^ny are the observed 
z/ fluxes. One can notice that, if the H transition is completely non adiabatic 
(Ph = 1), the NH and IH cases coincide. Thus, to see any effect due to the 
mass hierarchy, the H transition must be adiabatic, i.e. 9is has not to be too 
small. 

When we consider the effect of the Earth in the neutrino path to the detector, 
we must replace, in the detected ffux estimation in formulas (2-3-4-5), f/gj 
with Pjg (i = 1,2), the probability for the mass eigenstate z/j to be detected 
as fg (or Pg) after traveUing through the Earth [35], which depends on the 
solar oscillation parameters and on the travelled density profile through the 
Earth. We developed a complete 3-flavor calculation, describing the Earth 
interior as made of 12 equal density steps, following the Preliminary Reference 
Earth Model matter density proflle [36]. For each constant density step we 
compute the exact propagator of the evolution matrix and we get the global 



Pg=PHC/g^2^g°+(l-PHC/g'2)^° 



(2) 
(3) 



for normal hierarchy and 



F,^P^UlFl-r {l-P^Ul)Fl 



(4) 
(5) 
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Interactions with protons and iron nuclei Interactions with '^C nuclei 




Neutrino energy (MeV) Neutrino energy (MeV) 

Fig. 3. Theoretical cross sections of the neutrino CC interactions with protons and 
iron nuclei (left) and CC and NC interactions with carbon nuclei {right). 

amplitude matrix by multiplying the propagators of the traversed density 
layers, as described e. g. in [37]. 

More detail about the calculation of the probabilities Pie are given in Appendix 
A. 



4 LVD detector and the observable neutrino interactions 

The Large Volume Detector (LVD) in the INFN Gran Sasso National Labora- 
tory, Italy, consists of an array of 840 liquid scintillator counters, 1.5 m^ each, 
arranged in a compact and modular geometry; a detailed description can be 
found in [38,39]. The active scintillator mass is M = 1000 t. The counter are 
operated at energy threshold Sh — 5 MeV. To tag the delayed 7 pulse due to 
the neutron capture, all counters are equipped with an additional discrimina- 
tion channel, set at a lower threshold, ~ 1 MeV. The energy resolution is 
as/E = 0.07 + 0.23 • (E/MeV)-^-^ 



The observable neutrino reactions that occur in the scintillator are: 

(1) z/g p, e+n, (physical threshold Ep^ > 1.8 MeV) observed through a prompt 
signal from e+ above threshold £h (detectable energy E^ ~ £^p^ — 1.8 MeV 
+2 vfieC^), followed by the signal from the np,d ^ capture [E^ — 2.2 
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MeV) above £i and with a mean delay At ~ 185 /xs (as it comes out 
from estimation [40] and Montecarlo simulation [41]). The cross section 
for this reaction has been recalculated [42] with a better treatment of 
the 10 — 100 MeV region, i.e. the supernova neutrino energy. The cross 
section behavior with energy is shown in figure 3 {left plot, solid line). 
The total number of free protons in the scintillator is 9.34 ■ lO"'^. 

(2) z/e"'^^C,^^N e~, (physical threshold E^,^ > 17.3 MeV) observed through 
two signals: the prompt one due to the e~ above Sh (detectable energy 

~ E',^^ — 17.3 MeV) followed by the signal, above £h, from the /S"*" decay 
of ^^N (mean life time r = 15.9 ms). The efficiency for the detection of 
the beta decay product is 90% [43,44]. 

(3) z/e"'^^C,^^B e+, (physical threshold Ep^ > 14.4 MeV) observed through 
two signals: the prompt one due to the e"*" (detectable energy E^ ~ Ej^^ — 
14.4 MeV+2 mgC^) followed by the signal from the P~ decay of ^^B (mean 
fife time r = 29.4 ms). As for reaction (2), the second signal is detected 
above the threshold £h and the detection efficiency of the ^^B beta decay 
product is 75% [43,44]. 

(4) y« ^^C,u\ ^^C* {a = e,ii,T), (physical threshold E^ > 15.1 MeV) 
whose signature is the monochromatic photon from carbon de-excitation 
{E^ — 15.1 MeV), above £h, detected with 55% efficiency [45]. 

Cross sections for reactions (2), (3) and (4) are taken from [46] and shown 
in figure 3 (right). 

(5) V^cK e~, V^cK which yields a single signal, above £h, due to the recoil 
electron. Because of the low number of expected events (about a dozen) 
and the lack of a clear pattern for this interaction, we will not consider 
it in the following. 

The iron content in LVD (about 900 t) is concentrated in two components: 
the stainless steel tank (mean thickness: 0.4 cm) which contains the liquid 
scintillator and the iron module (mean thickness: 1.5 cm) which hosts a cluster 
of 8 tanks. Indeed, the higher energy part of the u flux can be detected also 
with the (iy) Fe interaction, which results in an electron (positron) that can 
exit iron and release energy in the scintillator. The considered reactions are: 

(6) i^e ^^Fe,^^Co* e~. The mass difference between the nuclei is Am„ = — 
mFe ^ 4 ^IpY ^Yie flrst Co aUowed state is at 1.72 MeV. Other 
allowed levels arc present in Cobalt, as shown in flg. B.l, whose energy 
i^ievei is 3.59, 7.2, 8.2, 10.59 MeV. Indeed, the electron kinetic energy 
is E^- = [Ey^ — — -Eievei — TTie) ■ Moreover, some gamma rays are 
produced in the interaction, depending on the excitation level considered. 

A full simulation of the LVD support structure and of the scintillator 
detectors has been developed in order to get the efficiency for electron and 
gammas, generated randomly in the iron structure, to reach the scintilla- 
tor with energy higher than S^. It is greater than 20% for E^, > 30 MeV 
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and grows up to 70% for E^, > 100 MeV. On average, the electron energy 
detectable is — 0.40 x E,y. The total number of iron nuclei is 9.22 • 10'^'^. 
(7) i/e^^Fe,^^Mn e^, the energy threshold is very similar to that of reaction 
(6). In this work, for simplicity, the same efficiency is assumed. 

The total cross section for reactions (6), (7) are taken respectively from [47] and 
[48] and plotted in figure 3 (/e/i), while the probability to select a particular 
Cobalt excitation level is taken from [49]. More detail about the neutrino-iron 
cross section, the Cobalt energy levels and the simulation of the interactions 
in the LVD detector are described in Appendix B. 

It is necessary to point out that, up to now, we calculated only v-Fc charged 
current interactions. The estimation of i/-Fe neutral current interaction cross 
section shows that they are roughly 30% of the CC ones [49] . They should be 
taken into account in future works. 

The number of all the possible targets present in the LVD detector is listed in 
table 1. 

Table 1 

Number of targets in the LVD detector. 



Target Type 


Contained in 


Mass (t) 


Number of targets 


Free protons 


Liquid Scintillator 


1000 


9.34 10^^ 


Electrons 


a 


1000 


3.47 10^2 


C Nuclei 


u 


1000 


4.23 10^^ 


Fe Nuclei 


Support Structure 


900 


9.71 10^° 



5 Expected neutrino signals 



The number of events detected during the supernova explosion is calculated 
as: 

oo 

TVe. = • / F{E;) ■ a{E,) ■ e{E,) dE, (6) 



where Nf is the number of target nuclei, F is the neutrino flux, a is the 
interaction cross section, E^, is the neutrino energy and e is the efficiency 
for the detection of the interaction products, where the effect of the detector 
energy threshold is included. 
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In the following wc show the effect of neutrino oscillations in the SN matter in 
the various interaction channels and the possible interplay among them. Then, 
the effect of the Earth matter is taken into account in the last subsection, 
considering only the inverse beta decay channel. 



5.1 Inverse beta decay 

The main interaction in LVD is the inverse beta decay (IBD) of electron 
antineutrinos. In figure 4 we show the energy spectra of the detected neutrinos 
in the case of no oscillation and in the case of adiabatic transition with NH 
and IH. We remind here that the non-adiabatic transition case (for both NH 
and IH) is coincident with the adiabatic NH case. 

In the case of oscillation, adiabatic, normal hierarchy, there is a contribution 
(sin^ 612) of the original higher-energy which gives rise to a higher average 
neutrino energy and, due to the cross section increase, to a larger number of 
detected events. The contribution is even higher (~ 1) if the transition is 
adiabatic and the hierarchy inverted, because the MSW resonance happens in 
the P sector. This results in a higher neutrino energy, as visible in figure 4, 
and in a larger number of events. This is clearly seen also in figure 5, where 
we show the number of z/g interactions with protons that can be detected in 
LVD as a function of the neutrino-sphere temperature. 



5.2 Charged current interactions with ^^C 

In figure 6 we show the expected number of (ug + i^e) charged current (CC) in- 
teractions with the ^^C nuclei. The two contributions have the same signature 
in the detector if one looks for two high threshold signals in a time window of 
about 100 ms, thus we consider them together. The conversion between the 
higher-energy non-electron neutrinos {u^, ) and the lower energy z/g, z/g , 
due to neutrino oscillation, increases the expected number of events. In the 
case of adiabatic transition the increase is even higher because at least one 
neutrino elicity state get a stronger contribution from the original i/^ (see eqq. 
2-5). 

A strategy to statistically determine the separate amount of Ue and Ue inter- 
actions, if a large number of CC interactions with ^^C is detected, is described 
in [50]. 



12 



5.3 Charged current interactions in the iron support structure 



An important contribution to the total number of events is also given by neu- 
trino interactions in the iron support structure of the LVD detector. Given the 
rather high effective threshold (about 10 MeV) and the increasing detection 
efficiency with the neutrino energy, they are concentrated in the high energy 
part of the spectrum [Ey > 20 MeV). Thus they are extremely sensitive to 
the neutrino energy spectrum and, indeed, to the oscillation parameters. 



Inverse beta decay interactions on protons 
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Fig. 4. Neutrino energy distribution in the interactions with p expected in LVD for 
three oscillation scenarios: no oscillation (solid line), adiabatic transition with NH 
(dashed), adiabatic transition with IH (dotted). The situation in the non adiabatic 
transition cases is identical to the adiabatic transition with NH case. The integral 
number of detected events is shown. 
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Fig. 5. Number of detectable Ve interactions with p expected in LVD as a function 
of the Vf. neutrino-sphere temperature. 



In figure 7 we show the dependence of the total number of detected (z/g+z/e) CC 
interactions with Fe to the i^e^sphere temperature , in the various oscillation 
scenarios. In figure 8 we show the contribution of [y^ + z/g) Fe interactions on 
the total number of events. For the chosen supernova and oscillation parame- 
ters they are about 17% of the total signal. Indeed, they have to be considered 
in an accurate estimation of the expected events. 
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5-4 Neutral current interactions on ^^C 



Neutral current interactions have the same cross section for all neutrino flavors, 
being thus insensitive to neutrino oscillations. Due to the high energy threshold 
(15 MeV) of the interaction, most of the detected events are given by the higher 
energy Ux- In principle, as shown in figure 9 (solid line), NC with ^^C could 
thus be used as a reference to identify the u^^sphere temperature. However, 
the expected number of events depends also on the value of other astrophysical 
parameters (as it will be discussed in section 6); for example, just changing 
the value of the other neutrino-sphere temperatures causes the variation in 
the number of events shown in figure 9, where two extreme values for the ratio 
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Fig. 6. Number of detectable {ug + Ug) CC interactions with ^^C expected in LVD 
as a function of the z^e neutrino-sphere temperature. 
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(i^e + k) CC interactions on Fe 
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Fig. 7. Number of detectable [y^ + fe) CC interactions with the iron of the support 
structure expected in LVD as a function of the v^. neutrino-sphere temperature. 

Ty^jTy^ are chosen: 1.1 (dotted) and 2. (dashed). 



5.5 Earth matter effect 



In order to measure the Earth matter effect at least two detectors in the world 
must detect the supernova neutrino signal, and one of them must be shielded 
by the Earth. In figure 10 the effect of Earth matter in the inverse beta decay 
interaction channel, that is the most abundant and the cleanest one, is shown. 
The nadir angle is Qn = 50°, which corresponds to neutrinos travelling only 
through the mantle. The Earth matter effect produces a decrease in the num- 
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Adiabatic Transition — Normal Hierarchy 
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Fig. 8. Neutrino energy distribution of the events occurring in the hquid scintillator 
(dashed), in the iron support structure (dotted) and their sum (solid) in the LVD 
detector. 

ber of detected neutrinos for particular neutrino energy, with a characteristic 
oscillating pattern. The effect is more relevant in the v than in the V channel, 

so the effect in reaction (1) is quite weak (the weakness of the effect also de- 
pends on the rather high Am^^^j = 8 • 10^^ eV^), but it could be detected if 
compared with a high statistic sample (i.e. with the Super-Kamiokande signal) 
or if a larger number of events is available (i.e. a closer supernova). 

Even if the Earth matter effect is very difficult to be detected, we want to 
mark its importance. The comparison between the distorted and undistorted 
spectrum cancels many uncertainties coming from the astrophysical parame- 
ters that, as described in detail in the next section, introduce a large variability 
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in the expected results. Then, the effect depends mainly on the 612 and Am^2 
oscillation parameters, which are now known with good precision. Moreover, 
even a null result carries many informations. For example, if the effect is not 
seen in the channel, it means that the transition is adiabatic and the hier- 
archy is inverted: in this case equation 5 becomes Fg = F^, and it does not 
depend on the transition probabilities Pie and P2e- An exhaustive discussion 
about the implications of the Earth matter effect in the supernova neutrino 
signal can be found e.g. in [51]. 




Fig. 9. Number of the detectable NC interactions of v and of all flavors with C 
expected in LVD as a function of the Px neutrino-sphere temperature. 
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Adiabatic Transition — Normal Hierarchy 
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Fig. 10. Effect of the Earth matter in the VeP, e^n signal in LVD, for neutrinos 
travelling through the Earth mantle (nadir angle 9n = 50°). 

6 Astrophysical parameter dependence 



As discussed in section 2, the emitted neutrino spectra depend on a large num- 
ber of astrophysical parameters, whose values are not uniquely established by 
the actual MonteCarlo calculations. Thus it is very important to take into ac- 
count the uncertainties in those parameters when estimating the total number 
of detected events. 



A summary of the values of the astrophysical parameters used in our calcu- 
lations is presented in table 2, together with an estimation of their range of 
variability, as attempted e.g. in [31]. The expected number of events in the 
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various LVD detection channels and the mean energy of the detected p 
events are shown in table 3. 



Astrophysical parameter 


Unit 


Chosen value 


Range of variability 


D: distance to the star 


kpc 


10 


0.2^20 


Ef,: total energy emitted in i/'s 


IQ53 g^rr 
o 


3. 


2. ^5. 


f,^^: fraction of E'b taken by Ug 




1/6 


1/10-1/4 


Tp^: Pe~sphere temperature 


MeV 


5. 


4.^7. 






0.8 


0.5^0.9 






1.5 


1.1^2. 


rj: pinching parameter 




0. 


0. H-2. 



Table 2 

Astrophysical parameters values used in the calculations and their assumed uncer- 
tainties. 





No Oscillation 


Non Adiabatic 


Adiabatic NH 


Adiabatic IH 


i^e P 


346. 


391. 


494. 


{Ep^) in Ue p 


25. MeV 


30. MeV 


37. MeV 


CC with 


8. 


22. 


29. 


27. 


CC with ^^Fe 


22. 


72. 


95. 


92. 


NC with 


27 



Table 3 

Expected results in the various LVD detection channels and in the mean energy of 
the detected z/g p events, calculated considering the chosen values of the astrophysical 
parameters, as given in table 2. 



Our approach will be the following: wc change the value of one parameter 
from the value listed in table 2 to the extreme lower and upper values, leaving 
the other astrophysical parameters unchanged and we show the fractional 
modification in the number of events and in the mean energy of the detected 
Pe in Pe P interactions. 

The distance to the supernova (D) and the total energy emitted as neutrinos 
(Eh) appear in the time-integrated neutrino spectrum (equation 1) as mul- 
tiplier factors {Nev oc Eb/D"^). Thus, their uncertainty can greatly modify 
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the expected number of events. However, the use of the ratios of appropriate 
observables (i.e. inverse beta decay vs neutral current with carbon) allow to 
factorize them away. The energy spectra is not distorted by a change in 
nor D, so < Ep^ > is unmodified. 



The energy taken by each neutrino flavor is usually considered as equipar- 
titioned; but differences up to a factor of two are allowed. If we consider 
fpe = Ue and the normahzation (/j,^ + fp^ + Af,yJ = 1, we can choose f^^ as 
the only independent parameter. We show in table 4 the fractional modifica- 
tions corresponding respectively to the values f^^ = 1/10 (— > fi,^ = 1/5) and 
U = 1/4 = 1/8). 





No Oscillation 


Non Adiabatic 


Adiabatic NH 


Adiabatic IH 


i^e P 


-40% , +50% 


-18% , +22% 


+20% , -25% 


(Ep^) in Pe p 


0% 


+7% , -6% 


0% 


CC with 


-40% , +51% 


+8% , -10% 


+13% , -16% 


+18% , -22% 


CC with s^Fe 


-40% , +50% 


+9% , -12% 


+ 14% , -17% 


+ 18% , -23% 


NC with 


+13% , -16% 



Table 4 

Fractional variations in the expected results if the energy fraction taken by is 
1/10 (left value) or 1/4 (right value), with respect to the chosen value 1/6. 



The values of the fg, i>e and neutrino-sphere temperatures determine the 

energies of the incoming neutrinos. They are parametrized in terms of the 
z/g-sphere temperature Tp^, and the ratios T^^/Tp^ and T^^/Tp^. The fractional 
variations when changing these three parameters are shown respectively in 
tables 5, 6, 7. 
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No Oscillation 


Non Adiabatic 


Adiabatic NH 


Adiabatic IH 




-20% , +35% 


-19% , +33% 


-17% , +29% 


{Ep^) in z/e p 


-19% , +37% 


-19% , +36% 


-19% , +37% 


CC with 


-56% , +164% 


-44% , +94% 


-43% , +85% 


-41% , +79% 


CC with 56Fe 


-55% , +193% 


-51% , +164% 


-51% , +161% 


-49% , +152% 


NC with 


-40% , +77% 



Tabic 5 

Fractional variations in the expected results if the i>(, neutrino-sphere temperature 
is equal to 4 MeV (left value) or 7 MeV (right value), with respect to the chosen 
value 5 MeV. 





No Oscillation 


Non Adiabatic 


Adiabatic NH 


Adiabatic IH 


l^e P 


0% 


(Ep^) in i/g p 


0% 


CC with ^-c 


-31% . +21% 


-5% . +3% 


0% . 0% 


- 1% . +2% 


CC with sepe 


-31% , +21% 


-3% , +2% 


0% , 0% 


-3% , +2% 


NC with 


-4% , +2% 



Table 6 

Fractional variations in the expected results if the ratio between the z^g and the 
neutrino-sphere temperature is equal to 0.5 (left value) or 0.9 (right value), with 
respect to the chosen value 0.8. 
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No Oscillation 


Non Adiabatic 


Adiabatic NH 


Adiabatic IH 




0% , 0% 


-9% , +10% 


-24% , +25% 


{Ep^) in z/e p 


0% , 0% 


-13% , +19% 


-25% , +31% 


CC with 


0% , 0% 


-43% , +54% 


-48% , +60% 


-50% , +59% 


CC with 56Fe 


0% , 0% 


-52% , +106% 


-57% , +116% 


-59% , +116% 


NC with 


-44% , +51% 



Tabic 7 

Fractional variations in the expected results if the ratio between the and the 
Ve neutrino-sphere temperature is equal to 1.1 (left value) or 2. (right value), with 
respect to the chosen value 1.5. 



The neutrino energy spectrum is presumably a black-body of the Fermi-Dirac 
type, but possible non-thermal effects are taken into account by introducing 
the parameter r] in equation 1, the so-called "pinching" factor. A distribution 
with 77 > and the same average energy is, in fact, suppressed at low and 
high energies. In table 8 we show the fractional differences in the results when 
considering t] = 1 01 t] = 2, m the simplified scenario where all the neutrino 
flavors are described by the same pinching factor. 





No Oscillation 


Non A(lial)atic 


A(lial)atic NH 


Adiabatic IH 


P 


+3% , +9% 


+3% , +9% 


+3% , +8% 


{EpJ in Ue P 


+2% , +5% 


+1% , +5% 


+2% , +5% 


CC with 


+7% , +21% 


+6% , +18% 


+6% , +18% 


+6% , +17% 


CC with ^'^^Fe 


+7% . +20% 


+6% . +19% 


+7% . +19% 


+6% . +18% 


NC with 


+6% , +17% 



Table 8 

Fractional variations in the expected results if the pinching parameter i] is equal to 
1 (left value) or 2 (right value), with respect to the chosen value r] = 0. 



We conclude that the sources of uncertainty in the astrophysical parameters 
which mostly affect the results are the partition of the available energy among 

the neutrino flavors and the values of the various neutrino-sphere tempera- 
tures. The largest variations in the expected signal (up to more than 100%) are 
hence due to the poor (and hard to get) theoretical knowledge of the physics of 
the gravitational collapse, which will be hopefully improved at the occurence 
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and detection of the next galactic supernova. With respect to the distance 
and the total released energy, which appear as a EiJ D"^ multiplier factor, 
their uncertainties certainly affect the signal, but, even if not constrained by 
the observation, they can be factorized away by using appropriate observables 
(e.g. the ratio of the NC events and the Pg P events). 



7 Summciry and conclusions 

The main aim of this paper was to show how neutrino oscillations affect the 
signal expected in the LVD detector at the occurence of the next galactic 
supernova. 

The LVD detector has been described in its main components. It is able to 
detect neutrinos of all flavors, by studying them in the various CC and NC 
channels. All the neutrino interactions that occur in the liquid scintillator as 
well as in the iron support structure have been studied in detail taking into 
account the neutrino energy threshold, cross section and detection efficiency. 

We assumed a galactic supernova explosion at a typical distance oi D — 
10 kpc, parametrized with a pure Fermi-Dirac energy spectrum (77 = 0) with 

a total energy i?5 = 3 ■ 10^^ erg and perfect energy equipartition /i,^ = /p^ = 
= 1/6; we fixed T,JT^^ ^ 1.5, TJT^^ = 0.8 and T^^ - 5 MeV. ' 

We considered neutrino oscillations in the standard three-flavor scenario. The 
MSW effect has been studied in detail for neutrinos travelling through the 
supernova matter. We also considered the distortion in the expected neutrino 
spectra induced by a possible path inside the Earth before their detection. 

For the chosen supernova parameters, it results that the expected number of 
events and their energy spectrum depend on the unknown oscillation param- 
eters: the mass hierarchy and the value of ^13. 

In particular, the inverse beta decay interactions (z/gj?, e'^n) are highly sensitive 
to the mass hierarchy: for adiabatic transition, the number of events increases 
of ~ 25% in the IH case, with respect to the NH one, since the detected 
completely come from the higher energy i/^- The mean energy of the detected 
positrons is correspondingly increased. 

The total number of (fg -|- Pg) CC interaction with ^^C nuclei is highly in- 
creased taking into account neutrino oscillations, because of their high energy 
threshold. For adiabatic transition the expected number of events is higher 
than the non adiabatic one, because at least one specie (between i/g or z/g) 
comes signiflcantly from the original and higher-energy Vy. in the star. How- 



24 



ever, if it is not possible to discriminate between and z/g, the normal and 
inverted hierarchy cases present similar results. Indeed, in the NH (IH) case, 
the increase in (i/e) is compensated by a decrease in {v^). 

The neutrino interactions with the iron of the support structure, which are 
studied in detail in this work, are also incread by the oscillations. The efficiency 
for the detection of the produced charged leptons and gammas in the active 
part of the detector has been calculated with a full simulation of the apparatus. 
The contribution of (z/g + Pe) Fe interactions can be as high as 17% of the total 
number of events (in the adiabatic NH case) and they contribute mostly to 
the high energy part of the spectrum. 

With respect to the previous detection channels, the number of NC interac- 
tions with ^^C nuclei does not depend on oscillations. In principle they could 
be used as a reference to identify the z/^^-sphere temperature. However, this is 
partly limited by the uncertainties in the other astrophysical parameters. 

We completed the calculations taking into account the effect of the passage of 
neutrinos through the Earth before their detection. This induces a character- 
istic modulation in the energy spectrum; however, given the expected number 
of events and the assumed oscillation parameters, the effect is quite weak. 

In conclusion, for the choice of the astrophysical parameters adopted in this 
work, the expected signal of neutrinos in the LVD detector from a supernova 
core collapse greatly benefits of the neutrino oscillation mechanism, practically 
in all the possible detection channels, especially if the transition is adiabatic 
and the hierarchy inverted (since in LVD the most relevant signal is given by 

However, being aware of the fact that the astrophysical parameters of the su- 
pernova mechanism are up to now not well defined, we performed the same 
calculations using different values of them. The resulting differences are in 
fact important; they are mainly due to the poor theoretical knowledge of the 
physics of the gravitational collapse. This will be hopefully improved after the 
occurence and detection of the next galactic supernova, to which the LVD ex- 
periment can give a significant contribution, thanks to its cabability to observe 
and measure neutrino events of several types. 
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Appendices 



A Oscillation probability in the Earth 



The neutrino flavor li/^) {a 
are related by: 



e,ii,T) and mass \uj) {i — 1,2,3) eigenstates 



(A.l) 



(for antineutrino U* should be replaced by U) where U is the 3x3 mixing 
matrix in vacuum 



Cl2 Ci3 



■S12 Ci3 



-Si2 C23 - C12 S23 -Sis e'^ C12 C23 - S12 S23 -Sis e*'' -§23 C13 

•S12 -S23 ~ C12 C23 Si3 e*"^ — C12 S23 ~ 512 C23 •S13 C23 Ci3 



(A.2) 



Since in the study of supernova neutrinos we are interested only in the v^. and 
survival probabilities, the angle ^23 and the CP-violating phase 5 do not 
matter. 

The flavor eigenstates evolution in matter is governed by the following equa- 
tion: 



.0? 1 
dt ' ' ^2E 



(A.3) 



where is the diagonal matrix of the squared neutrino masses and Vcc is 
the matter induced potential 3x3 matrix 



^1 0^ 








(A.4) 



with Gp the Fermi constant and Ng the electron density in the crossed matter. 
Diagonalizing H we get 

H = UmDUl (A.5) 
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where D is the diagonal matrix of the cigcnvahics in matter and Um is the 
orthogonal matrix with the mass eigcnstatcs in matter as column. We thus 
define jt'™), the neutrino eigenstates for the propagation in matter of constant 
density N^, as 

WD = (Ul),^ Wa) (A.6) 



and the evolution equation is 



The mass eigenstates at the time t become 

krW)=e-^* kr(0)) (A.8) 

and the flavor eigenstates at the time t are 

Wait)) = UrnB-'^^'Ul |z/„(0)) = S{t,0) MO)) (A.9) 



where S{t,0) = Ume f/^ is the propagator of the flavor eigenstate 
from the time t = to the time t, in the matter of constant density Ng. 

This result is valid if the neutrinos travel through a single layer of constant 
density described by Ne- Describing the density distribution as a series of 
steps, each of constant density, as it is for the Earth interior, we only need to 
replace the propagator S{t, 0) with the product of many propagators, one for 
each crossed density layer. Suppose, for example, that the neutrinos enter the 
layer 1 with density at the time ti and exit at t2 going through the layer 2 
with density until the time ^3 when they are detected. In this case we can 
decompose the propagator as 

^(^3, h) = ,^2(^3, ^2) • Si{t2, h) (A.io) 



where each propagator has to be calculated as in eq. A.9 considering the 
particular density of the corresponding layer. 

In our calculations the Earth interior has been divided into 12 layers of con- 
stant density, following the PREM model [36]; they are shown in flgure A.l. 
Two main regions can be deflned: the Core (radius < 3500 km, nadir angle 
On < 33°) and the Mantle (outside). 
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Fig. A.l. Density of each layer inside the Earth in term of their radius (lower X 
axis). The nadir angle (0^) corresponding to each shell radius is shown in the top 
level of the figure. The vertical dashed line at 3480 km {6n = 33°) represents the 
division between the Core and the Mantle inside the Earth, where a large change 
in density occurs. 

Now, since we are interested in calculating the probability for the mass eigen- 
state j, coming from the supernova, to be detected as flavor (3, we get 

P(j^P) - \{Mo)Wjm' - \ E ul s{t,o)f,^\' (A.ii) 

a=e,iJ,,T 



The results of the calculation of are shown in figure A. 2. In the left plot 
the nadir angle is 20° and the neutrinos go through both the Mantle and the 
Core. The interference between the two main density layers give rise to the 
complicated behavior of the probability Pig. In the right plot 6n = 50° and 
the path is in the mantle only. The traversed density is thus almost constant 
and Pie becomes more regular. An animated version of figure A. 2, showing 
continuously what happen when the nadir angle changes, can be found in 
[52]. 

Similar plots can be obtained for the probability. However, since P^^ is very 
close to zero, for our purposes P2e = 1 — Pie is a very good approximation. 
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Earth matter effect — nadir angle = 20° 



Earth matter effect — nadir angle = 50" 
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Fig. A. 2. Pie is the probability that the neutrino mass eigenstate z^i, coming from 
the supernova, is detected as I'e after travelhng through the Earth, as a function of 
the neutrino energy. The red sohd hne is for antineutrinos, while the blue dotted 
one is for neutrinos. The horizontal black solid line represents the value of Pie in the 
case of no Earth matter effect: Pie = cos^(0i2). In the left plot the case of neutrinos 
travelling through the Mantle and the Core is shown (0„ = 20°), while in the right 
plot 6n = 50° and the neutrinos go only through the Mantle, where the density is 
almost constant. 

When the neutrino path is only through the Mantle, the crossed density can 
be approximated as a single constant density step (p ~ 4 g/cm^). In this 
particular case the Pie probabilities are exactly described by [31]: 



Pie = 



Ae sin^ 9i2 cos^ 9 



12 



;i + -4ecos2 Ou 



sm 



1 + -4£cos2^i2 (A.12) 



where 

^ V2GFNe 

^ Aml2/2E 
For Ue, just replace Ou - 



90° - 6i2. 



p/{A g/cm^) • re/(0.5) ■ E/(20 MeV) 



Amf2/(8 ■ 10-5 eV^ 
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B Neutrino-iron interactions in LVD 



We consider the reaction ^^Fe,^^Co* e~ . It leads to excitation of analog 0+ 
and Gamow-Teller 1+ giant resonances (AR and GTR, respectively) in ^^Co 
nucleus. The excitation of the AR is connected with Fermi transition. The 
ground state quantum numbers of ^^Co are 4+, therefore the corresponding 
cross section with ^^Co^.s. in the final state is small, compared to AR and 
GTR excitation. 

A simplified Cobalt nucleus level structiue is shown in figure B.l, where the 
O"*" AR is 3.59 MeV over the ground state and the GTR take the other energy 
levels. Therefore, in our simulation, we take into account 5 possible channels 
[49]: 

• AR) Cobalt is excited to the 0+ AR, at 3.59 MeV. The products of the 

interaction are: 

• the electron, with kinetic energy E^,- = Ey^ — ( A^^ + i^ievei + i^e) = — 
8.156 MeV, where ^levei = 3.59 MeV and A„„ = m^°-mf = 4.055 MeV, 

• one 1.87 MeV gamma, 

• a gamma cascade, whose total energy is 1.72 MeV 

• GTl) Cobalt is excited to the first Gamow-Teller resonance, at 1.72 MeV. 
The products of the interaction are: 

• the electron, with kinetic energy E^- — E^^ — 6.286 MeV, 

• a gamma cascade, whose total energy is 1.72 MeV 

• GT2) Cobalt is excited to the second GTR, at 7.2 MeV. The products of 
the interaction are: 

• the electron, with kinetic energy E^- — E^,^ — 11.766 MeV, 

• one 3.61 MeV gamma, 

■ one 1.87 MeV gamma, 

• a gamma cascade, whose total energy is 1.72 MeV 

• GT3) Cobalt is excited to the third GTR, at 8.2 MeV. The products of the 
interaction are: 

• the electron, with kinetic energy E^- — E^,^ — 12.766 MeV, 

■ one 4.61 MeV gamma, 

• one 1.87 MeV gamma, 

• a gamma cascade, whose total energy is 1.72 MeV 

• GT4) Cobalt is excited to the fourth GTR state, at 10.6 MeV. The products 

of the interaction are: 

• the electron, with kinetic energy E^.- = E^^ — 15.166 MeV, 

• one 7. MeV gamma. 
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one 1.87 MeV gamma, 

a gamma cascade, whose total energy is 1.72 MeV 
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Fig. B.l. Excitation levels of the 56 Cobalt nucleus. 

The neutrino- iron partial cross sections are shown in table B.l for each con- 
sidered excitation channel; they are taken from 



For each neutrino energy the excitation channel is sampled accordingly to its 
relative weight to the total cross section. The reaction products are generated 
uniformly in the LVD iron support structure. The electron and the gammas 
are simulated starting from the same generation point and their directions are 
chosen uniformly in the whole solid angle, without any correlation between 
them. 



We used a GEANT3 simulation of the LVD detector, where the liquid scintil- 
lator and the iron support structure are described in detail. The particles of 
the electromagnetic showers are tracked through the various materials until 
their energy is smaller than 100 keV, the lowest possible value allowed by the 
simulation program. 

We define the efficiency as the ratio of the number of events where at least one 
scintillation counter detects a signal over its threshold and the total number of 
generated events. The energy resolution of the detector is taken into account 
in the simulation. Setting the energy threshold of all the scintillation counters 
to 5 MeV, the resulting efficiency is shown in figure B.2. 
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Table B.l 

Neutrino-iron cross section for the various excitation level considered. 



Neutrino energy (MeV) 


Cross section (10 ^° cm^) 




GTl 


AR 


GT2 


GT3 


GT4 


10 


8.83 10-3 


7.29 10-3 


- 


- 


- 


20 


9.58 10-2 


1.96 10-1 


1.23 10-2 


5.81 10-2 


6.76 10-2 


30 


2.72 10-1 


6.26 10-1 


5.49 10-2 


2.96 IQ-i 


5.49 10-1 


40 


5.36 10-1 


1.29 10° 


1.27 10-1 


7.10 10-1 


1.48 10° 


50 


8.86 10-1 


2.19 10° 


2.28 10-1 


1.30 10° 


2.84 10° 


60 


1.32 10° 


3.32 10° 


3.57 10-1 


2.06 10° 


4.63 10° 


70 


1.84 10° 


4.68 10° 


5.15 10-1 


2.99 10° 


6.85 10° 


80 


2.45 10° 


6.26 10° 


7.01 10-1 


4.09 10° 


9.50 10° 


90 


3.14 10° 


8.07 10° 


9.15 10-1 


5.36 10° 


1.26 10+1 


100 


3.90 10° 


1.01 lOi 


1.16 10° 


6.79 10° 


1.60 10+1 



The total energy detected in the liquid scintillator is very weakly correlated 
to the neutrino energy, as shown in the scatter plot of figure B.3. On average, 
the total energy detectable is E'^ ~ 0.4 x E^^,, but the spread over the mean 
value is very large. 
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Fig. B.2. Detection efficiency of neutrino-iron interaction for an energy threshold of 
the scintillation counters of 5 MeV. 
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